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ABSTRACT 

We provide evidence that the obliquities of stars with close-in giant planets were initially nearly random, and 
that the low obliquities that are often observed are a consequence of star-planet tidal interactions. The evidence 
is based on 14 new measurements of the Rossiter-McLaughlin effect (for the systems HAT-P-6, HAT-P-7, HAT- 
P-16, HAT-P-24, HAT-P-32, HAT-P-34, WASP-12, WASP-16, WASP-18, WASP-19, WASP-26, WASP-31, 
G1436, and Kepler-8), as well as a critical review of previous observations. The low-obliquity (well-aligned) 
systems are those for which the expected tidal timescale is short, and likewise the high-obliquity (misaligned 
and retrograde) systems are those for which the expected timescale is long. At face value, this finding indicates 
that the origin of hot Jupiters involves dynamical interactions like planet-planet interactions or the Kozai effect 
that tilt their orbits, rather than inspiraling due to interaction with a protoplanetary disk. We discuss the status 
of this hypothesis and the observations that are needed for a more definitive conclusion. 

Subject headings: techniques: spectroscopic — stars: rotation — planetary systems — planets and satellites: 
formation — planet-star interactions 



1. INTRODUCTION 

Exoplanetary science has been full of surprises. One of the 
biggest surprises emerged at the dawn of this field: the exis- 
tence of "hot Jupiters" having orbital distances much smaller 
than an astronomical unit (AU). It is thought that giant planets 
can only form at distances of several AU from their host stars, 
where the environment is cooler and solid particles are more 
abundant, facilitating the growth of rocky cores that can then 
attract gaseous envelopes from the protoplanetary disk. 

Different mechanisms have been proposed which might 
transport giant planets from their presumed birthplaces in- 
ward to where we find them. Among the differences between 
the proposed mechanisms is that some of them would alter 
the planet's orbital orientation and thereby change the rela- 
tive orientation between the stellar and orbital spin (Naga- 
sawa et al. 2008; Fabrycky & Tremaine 2007; Naoz et al. 
2011), while others would preserve the relative orientation 
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(Lin et al. 1996), or even reduce any primordial misalignment 
(Cresswell et al. 2007). For this reason, measuring the stel- 
lar obliquity — the angle between stellar and orbital axes — has 
attracted attention as a possible means of distinguishing be- 
tween different theories for the origin of hot Jupiters. 

The stellar obliquity is an elusive parameter because the 
stellar surface needs to be at least partially resolved by the ob- 
server. If the system exhibits eclipses or transits, then it is pos- 
sible to detect anomalies in the spectral absorption lines of the 
eclipsed star, which have their origin in the partial blockage of 
the rotating photosphere. The precise manifestation of the ro- 
tation anomaly depends on the angle between the projections 
of the stellar rotation axis and orbital axis of the occulting 
companion. Credit for the first definitive measurements has 
been apportioned to Rossiter (1924) and McLaughlin (1924), 
after whom the effect is now named. Queloz et al. (2000) were 
the first to measure the Rossiter-McLaughlin (RM) effect for a 
planet-hosting star, finding a low obliquity. Since then many 
systems have been studied, including those with misaligned 
stars (Hebrard et al. 2008), retrograde orbits (Winn et al. 2009; 
Triaud et al. 2010), a hot Neptune (Winn et al. 2010c; Hi- 
rano et al. 201 lb) and even a circumbinary planet (Winn et al. 
2011b). 

Winn et al. (2010a) found a possible pattern in the hot- 
Jupiter data, namely, host stars hotter than r e g w 6250 K tend 
to have high obliquities, while cooler stars have low obliqui- 
ties. Schlaufman (2010), using a different method, also found 
a preponderance of high obliquities among hot stars. Winn 
et al. (2010a) speculated that this pattern was due to tidal in- 
teractions. Specifically, it was hypothesized that hot Jupiters 
are transported inward by processes that perturb orbital incli- 
nations and lead to a very broad range of obliquities. Cool 
stars ultimately come into alignment with the orbits because 
they have higher rates of tidal dissipation, due to their thick 
convective envelopes. Hot stars, in contrast, lack thick con- 
vective envelopes and are unable to reorient completely on 
Gyr timescales. 
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If this interpretation is correct, then not only should fur- 
ther measurements be consistent with the hot/cool pattern, 
but also the degree of alignment should be found to correlate 
with the orbital period and planet-to-star mass ratio, param- 
eters which also strongly affect the rate of tidal dissipation. 
Since the study by Winn et al. (2010a) the number of sys- 
tems with RM measurements has nearly doubled. Here we 
present the results of an additional 14 observations, as well as 
a critical review of other published measurements (including 
thorough re-analyses in three cases). We can now attempt a 
comparison between the measured obliquities and the theoret- 
ical distribution of obliquities that one would expect if tides 
were an important factor. We also refer the reader to Hansen 
(2012), who recently performed a comparison with a similar 
motivation, without the benefit of the new RM measurements 
presented here, but using more sophisticated theoretical mod- 
els and drawing qualitatively similar conclusions. 

This paper has two main parts. The first part is observa- 
tional. We describe our new observations of the RM effect 
(Section 2), our analysis method (Section 3), and the details 
of individual systems analyzed here (Section 4). The second 
part (Section 5) considers the distribution of obliquities, seeks 
evidence for the expected signatures of tidal effects, and con- 
siders the implications for the origin of hot Jupiters. We sum- 
marize our results in Section 6. 

2. OBSERVATIONS 

The 14 new observations presented in this paper were con- 
ducted with the Keck I telescope and its High Resolution 
Spectrograph (HIRES; Vogt et al. 1994), for northern objects, 
and with the Magellan II (Clay) 6.5 m telescope and the Planet 
Finder Spectrograph (PFS; Crane et al. 2010) for southern ob- 
jects. Table 1 is a log of the observations. To derive the rela- 
tive RVs we compared the spectra observed through the iodine 
cell with the stellar template spectrum multiplied by an iodine 
template spectrum. The velocity shift of the stellar template 
as well as the parameters of the point-spread function (PSF) of 
the spectrograph are free parameters in this comparison. The 
velocity shift of the template that gives the best fit to an ob- 
served spectrum represents the measured relative RV for that 
observation. In particular we used codes based on that of But- 
ler et al. (1996). The RVs used in this paper for all system are 
presented in the Table 5. 

3. ANALYSIS OF THE ROSSITER-MCLAUGHLIN EFFECT 

In the following sections we describe our strategy to ob- 
tain projected obliquities from the observed radial velocities 
(RVs). Our goal was to perform a relatively simple and homo- 
geneous analysis and thereby facilitate comparisons between 
systems in the sample. For some systems for which we ob- 
tained new data, RM measurements have already been pub- 
lished by other authors. In order to obtain independent mea- 
surements and compare the results, in most cases we did not 
include the previous data in our analysis. 

Observing the Rossiter-McLaughlin (RM) effect allows one 
to measure the projected angle between the orbital and stel- 
lar spins (A) and the projected stellar rotation speed (vsini*)- 
Here v indicates the stellar rotation speed and z* the inclination 
of the stellar spin axis towards the observer. Observing the 
RM effect provides information on the transit geometry be- 
cause the distortion of the stellar absorption line — specifically 
its net redshift or blueshift — depends on the RV of the hidden 
portion of the stellar photosphere. 

For example if the obliquity were low, then the planet would 



TABLE 1 
Observation Loo 



System 


Observation night 
UT 


Spectrograph 


HAT-P-6 


5/6 November 2010 


HIRES 


TTAT n T 

HA1-P-7 


23/24 July 2010 


HIRES 


TUT n 1 £. 

HA1-P-16 


26/2/ December 2010 


HIRES 


T I AT" T> H A 

HAI-P-24 


2//28 September 2010 


HIRES 


HAT-P-32 


5/6 December 2011 


HIRES 


HAT-P-34 


2/3 September 2011 


HIRES 


WASP- 12 


1/2 January 2012 


HIRES 


WASP- 16 


3/4 July 2010 


PFS 


WASP- 18 


8/9 October 201 1 


PFS 


WASP- 19 


20/21 May 2010 


PFS 


WASP-26 


17/18 August 2010 


PFS 


WASP-31 


12/13 March 2012 


HIRES 


G1436 


24/25 April 2010 


HIRES 


Kepler 8 


7/8 August 2011 


HIRES 



begin the transit on the approaching half of the photosphere. 
The blockage causes the absorption lines to appear slightly 
redshifted. During the second half of the transit, the reverse 
would be true: the anomalous RV would be a blueshift. In 
contrast, if the planet's trajectory were entirely over the red- 
shifted half of the star, then the RM effect would be a blueshift 
throughout the transit. For a retrograde configuration, the RM 
effect would first lead to a blueshift and then a redshift. 

In sections 3.1-3.5 we discuss some specific effects which 
influence the RM signal, and how we included these in our 
modeling. Section 3.6 gives the details of our quantitative 
model of the RM effect, and the prior information about each 
system that is used in the analysis. 

3.1. Line broadening 

In addition to stellar rotation, microturbulence and macro- 
turbulence will affect stellar absorption lines and should be 
considered in the description of the RM effect. Microturbu- 
lence is well described by a convolution of the rotationally- 
broadened line with a Gaussian function. Describing macro- 
turbulence is more complex, as it depends on the angle be- 
tween the line of sight and the local surface normal. Near the 
center of the stellar disk, the Doppler shifts are produced by 
vertical ("radial") motions, while near the limb, the Doppler 
shifts are produced by motions along the stellar surface ("tan- 
gential" motions). The net effect is obtained by spatial in- 
tegration over the visible hemisphere. Hirano et al. (2011c) 
provided a semi-analytic description of this effect, which we 
used in our model. The macro turbulent field is parameterized 
by £rt, the average value of the radial and tangential velocities 
(assumed to be equal). 

What value of £rt should be used to model the RM effect 
for a given system? Measuring macroturbulence from the line 
profile is challenging, as its measurement is correlated with 
the measurement of vsin/ + . In principle the RM effect it- 
self could be used to measure macroturbulence. However for 
most of the systems analyzed here the data sets are not precise 
enough to allow £rt to be a free parameter. If the macrotur- 
bulence cannot be determined by the RM effect, one might 
think its value also does not matter for the analysis. How- 
ever, as described below, we used prior information on v sin 
as a constraint on our model, and this information can only 
be properly incorporated if macroturbulence is taken into ac- 
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Fig. 1, — Illustration of the influence of rotation, PSF broadening and macroturbulence, differential rotation, and convective blueshift on the RM 
effect, (a) The effect of rotation. The chosen parameters are A = 40°, vsini* = 3 km s -1 , Rp/i?* = 0.12 and an impact parameter of 0.2. The x-coordinate is the 
planet's location relative to its position at inferior conjunction, (b) The effects of instrumental broadening (2.2 kms~') and macroturbulence (£rt = 3 kms~'). (c) 
The effect of solar-like differential rotation, (d) The effect of a solar-like convective blueshift. (e-h) The corresponding models for the stellar absorption lines. 
The red region correspond to the light lost during the particular transit phase indicated by a red dot in the upper panel. For the lower panels the radius of the 
planet was doubled, for improved visibility of the missing velocity components, (i) The combined model including all aforementioned effects. The gray line is 
the model from panel (a). 



count. We therefore made use of the relation given by (Gray 
1984, eq. 4), 

Crt = 3.95 7k- 19.25, (1) 

where £rt is in km s" 1 and r en - in thousands of Kelvin. Gray 
(1984) gives an uncertainty of 0.2 kms" 1 for this relationship 
between T e ff and £ RT . We adopt a more conservative uncer- 
tainty estimate of 1.5 kms" 1 . We decided to use the rela- 
tionship by Gray (1984) for macroturbulence instead of the 
relationship by Valenti & Fischer (2005), because the latter 
authors derived the relationship between T e s and (rt for F 
stars by extrapolation from later spectral types, while Gray 
(1984) obtained macroturbulence values for hotter stars di- 
rectly. Many of the stars in our sample have a r et f greater then 
6000 K. 

Additional contributions to the linewidth come from col- 
lisional broadening, and Zeeman splitting, which are small 
compared to the effects described above. We account 
for them, at least approximately, by convolving the disk- 
integrated line profiles with a Lorentzian function of width 
1 kms" 1 . 

3.2. Convective blueshift 

Upward-flowing material in a convective cell is hotter and 
more luminous than downward-flowing material. This leads 
to a net blueshift in the disk-integrated light, known as the 
convective blueshift (CB). The CB is strongest in light re- 
ceived from the center of the disk, and weaker near the stellar 
limb. Disk integration leads to an overall Doppler shift and an 
asymmetry in the stellar absorption lines. The overall Doppler 
shift is of order 1 km s" 1 , but this is unimportant for our pur- 
poses, as we are only interested in relative RVs. However, 
the motion of the transiting planet will cause time variations 
in the net CB of the exposed surface of the star. The time- 
varying component of the CB effect is of order 1 m s" 1 for 



transits of late-type stars (Figure 1 ). This is a small effect but 
for completeness it was included in our model. 

To describe the effect of the CB, we used a numerical model 
based on work by Shporer & Brown (2011). For the stars in 
our sample with T e ff < 6000 K we assumed the net convective 
blueshift to be similar to that of the Sun: an effective radial 
velocity of the photosphere of 500 m s" 1 . For hotter stars we 
assumed an effective radial velocity of 1000 ms" 1 . The sin- 
gle exception to this rule was Gl 436, which is much cooler 
than the other stars (r e ff = 3350 ± 300 K). For this system we 
assumed a smaller blueshift of 200 m s" 1 . In principle one 
could estimate the overall CB for each system based on the 
observed line bisectors, but we did not pursue this approach 
because our first-order analysis revealed that the CB effect is 
relatively unimportant. 

3.3. Differential rotation 

For this study we neglected the possible influence of differ- 
ential surface rotation on the RM effect. Albrecht et al. (2012) 
concluded that modeling differential rotation is justified only 
if the transit chord spans a wide range of stellar latitudes, the 
RM effect is detected with a very high signal-to-noise ratio, 
and excellent prior constraints are available for the limb dark- 
ening (not only as it affects the continuum but also the absorp- 
tion lines). These conditions are not met simultaneously for 
any of the systems in our sample. See Albrecht et al. (2012) 
for further discussion and an application of a model with dif- 
ferential rotation to the WASP-7 system. 

3.4. RV measurements during transits 

A model which aims to simulate RV measurements made 
during transits needs to take into account that the measured 
RVs represent the output of a Doppler-measuring code. These 
codes are complex but in essence they locate the peak of 
a cross-correlation between a stellar template spectrum ob- 
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TABLE 2 
Prior knowledge 



System 


Jitter 


T c - 2 400 000 


r 41 


T-n 


p/ * 


u\ + u.2 


Crt 


v sin i+ 


ret. 




(ms -1 ) 


(B JDTrm ) 


(hr) 


(hr) 






(kms 4 ) 


(kms~') 




HAT-P-6 


7.0 


54035.67652±0.00196 


3.506±0.041 


0.45 1 ±0.053 


0.0934±0.0005 


0.67±0.10 


5.52±1.50 


8.7±1.5 


1 


HAT-P-7 


6.0 


55401.91904±0.000()2 


3.941 ±0.002 


0.373±0.004 


0.0777±0.0001 


0.70±0.10 


5.18±1.50 


3.8±1.5 


2 


HAT-P-16 


4.0 


55027.59293±0.OO574 


3.062±0.031 


0.360±0.067 


0.1071±0.0014 


0.71±0.10 


4.88±1.50 


3.5±1.5 


3 


HAT-P-24 


12.0 


55216.97667±0.00071 


3.694±0.019 


0.338±0.029 


0.0970±0.0012 


0.70±0.10 


5.22±1.50 


10.0±1.5 


4 


HAT-P-32 


42.0 


54420.44637±0.00069 


3.108±0.007 


0.413±0.010 


0.1508±0.0004 


0.69±0.10 


4.96±1.50 


20.7±1.5 


5 


HAT-P-34 


25.0 


55431.59629±0.00126 


3.492±0.038 


0.290±0.062 


0.0801±0.0026 


0.69±0.10 


5.32±1.50 


24.0±1.5 


6 


WASP- 12 


4.5 


54508.97605±0.00146 


3.001 ±0.037 


0.324±0.040 


0.1119±0.0020 


0.71±0.10 


5.10±1.50 


2.2±1.5 


7.8 


WASP- 16 


10.0 


54584.42952±0.0U509 


1.934±0.031 


0.504±0.096 


0.10/9±0.0012 


0. 74±0.10 


A 1 O _L 1 C A 

4. 18±1.50 


3.0±1.5 


9 


WASP- 18 


8.0 


54221.48163±0.00155 


2.191±0.012 


0.218±0.034 


0.0970±0.0010 


0.68±0.10 


5.26±1.50 


11.0±1.5 


10 


WASP- 19 


20.0 


55168.96879±0.00011 


l.572±0.007 


0.324±0.024 


0.1425±0.0014 


0.76±0.10 


3.87±1.50 


4.0±2.0 


11,12 


WASP-26 


10.0 


55228.38916±0.00075 


2.383±0.043 


0.057±0.125 


0.101 1±0.0017 


0.74±0.10 


4.56±1.50 


2.4±1.5 


13,14 


WASP-31 


6.0 


55209.71890±0.00280 


0.110±0.00l 


0.029±0.004 


0.1 271 ±0.00 11 


0.69±0.10 


5.11±1.50 


7.9±1.5 


15 


G1436 


1.5 


54235.83624±0.01221 


1.0l2±0.010 


0.213±0.028 


0.0825±0.0078 


0.87±0.10 


0.56±1.50 


+ 5 


this work 


Kepler-8 


10.0 


55781.9061 1±0.00032 


3.257±0.010 


0.557±0.029 


0.0948±0.0006 


0.71±0.10 


4.97±1.50 


10.5±1.5 


16 


HAT-P-2 


0.0 


54387.49375±0.00075 


4.289±0.031 


0.338±0.072 


0.0723±0.0006 


0.70±0.10 


5.09±1.50 


20.8±1.5 


17 


HD 149026 


2.5 


54456.78835±0.00080 


3.230±0.150 


0.227±0.024 


0.0507±0.0009 


0.72±0.10 


4.89±1.50 


6.0±1.5 


18,19 


HD 209458 


2.0 


51659.93742±0.00002 


3.072±0.003 


0.438±O.0O8 


0.1211±0.0001 


0.71±0.10 


4.64±1.50 


4.5±1.5 


20 
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(2012); (7) Maciejewski et al. (201 1); (8) Hebb et al. (2009); (9)Lister et al. (2009); (10) Hellier et al. (2009); (11) Hebb et al. (2010); (12) Hellier et al. (201 1); 
(13) Smalley et al. (2010); (14) Anderson et al. (2011a); (15) Anderson et al. (2011c); (16) Jenkins et al. (2010); (17) Pal et al. (2010); (18) Carter et al. (2009); 
(19) Sato et al. (2005); (20) Laughlin et al. (2005) 



tained outside transit, and a spectrum observed during transit. 
Depending on the system parameters these RVs can signifi- 
cantly differ from RVs representing the first moment of the 
absorption line. The later are often used in RM work poten- 
tially leading to systematic errors in the derived parameters 
(Hirano et al. 2010, 2011c). 

Hirano et al. (2011c) presented a semi-analytic description 
of the shift in the cross-correlation peak as a function of the 
transit parameters, including the dependence of the RVs on 
the stellar rotation velocity and obliquity, the microturbulent 
and macroturbulent velocities, the differential rotation profile, 
and the PSF of the spectrograph. We used this model, after ex- 
tending it to include the convective blueshift. We also tested 
the results of the code by Hirano et al. (201 lc) with the fully 
numerical disk-integration code of Albrecht et al. (2007), in- 
cluding the same physical effects, finding good agreement. 
Figure I shows the influence of the effects discussed in this 
and the previous sections on the stellar absorption lines. 

3.5. Other RV variation sources 

In addition to the RM effect we must also model the 
changes in radial velocity due to the star's orbital motion. 
Over the course of the transit night, this motion can usually 
be represented as a linear function of time, and parameterized 
by a constant acceleration. For ease of comparison with ex- 
isting orbital solutions, the out-of-transit velocity trend may 
also be parameterized by the orbital velocity semi-amplitude 
(K+), for fixed values of the other orbital parameters (period 
P, eccentricity e, argument of pericenter lu). 

The question arises whether to allow to be a completely 
free parameter in our analysis, or whether to use a prior con- 
straint based on a previously published orbital solution. The 
answer is not obvious because stars exhibit intrinsic RV noise 
that may have different amplitudes on different timescales, 
ranging from a few minutes to days. Short-timescale noise 
will simply degrade our measurement accuracy, but noise on 



timescales longer than ^6 hr (such as the noise produced by 
rotating starspots) will introduce trends in the RV signal over 
the course of the transit night. Without a specific model for 
the frequency content of the intrinsic RV noise, it is difficult to 
combine the data obtained sporadically over days or months 
with the data obtained with much finer time sampling on a 
single night. Specifically, if we would choose to constrain 
K+ (and therefore the transit-night velocity gradient) based on 
the orbital solution, there is a risk of introducing a bias in our 
results because of actual transit-night velocity gradient could 
be affected by starspot variability or other intrinsic RV noise 
with timescales longer than a few hours (Bouchy et al. 2008; 
Albrecht et al. 2011b, 2012). 

One might proceed by downweighting the constraint on K+ 
by some amount deemed to be appropriate, or fitting all of 
the RV data using a frequency-dependent noise model. For 
this work we placed a high premium on simplicity and homo- 
geneity, and decided to decouple the transit-night data from 
the rest of the data. We allowed K+ to be a completely free pa- 
rameter and determine the velocity gradient from the transit- 
night data alone. This was done for all the systems in our 
sample, even those where we find no obvious sign of intrin- 
sic RV noise. This leads to larger uncertainties, and for some 
systems we find wider confidence intervals for A than other 
researchers using similar data or even noisier data. While in 
some systems this might be an unnecessary precaution, the 
stars' characteristics are in general not well enough known to 
make a principled case-by-case decision. Table 3 gives the 
values of K* from published orbital solutions, along with the 
the values of K+ derived from our model which was fitted to 
the transit-night data. 

For those readers who do not agree with our reasoning on 
this point or who are simply curious about the effect of im- 
posing a prior on K+, Table 3 and Figure 2 also provide the 
results for A and vsin/* obtained by imposing a prior on K+ 
based on the orbital solution. 
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TABLE 3 
Results for rotation parameters 



Final results With prior on K^' 



Oy alClll 


* /i) cin i cin A 
A/ 1 ftlll t-j, Sill /\ 


Ay I alll 1 * LAJft /\ 




\ 


I\ . 11 ill 1M I 


K 1 1 1^1"^ tl ItY 3 
I\ * IHL1 JIU1 L 




\ 


Ref.' 1 




(kins -1 ) 


(kms -1 ) 


(km s _1 ) 


c c '^ 


(kms _l ) 


(km s _1 ) 


(km s _1 ) 


V ) 




HAT-P-6 


0.702=1=0.269 


-2.689=1=0. 131 


7.8±0.6 


165±6 


188±20 


115.5=b4.2 


7.0=1=0.6 


175±4 


1 


HAT-P-7 


0.695±0.361 


-1.476±0.196 


2.7=b0.5 


155±37 c 


214±5 


213.5=bl.9 


2.7=1=0.4 


155±14 


2 


HAT-P-16 


-0.051±1.413 


1.327±0.364 


3.1=1=1.0 


9+55 
-46 


536=b60 


53I.1±2.8 


2.7=1=0.8 


-6±37 


3 


HAT-P-24 


1.158=1=0.851 


3.087±0.297 


11.2=1=0.9 


20=1=16 


129±41 


83±3.4 


ll.OiO.9 


14±16 


4 


HAT-P-32 


4.517=1=0.168 


0.396=1=0.117 


20.6=1=1.5 


85.0=1=1.5 


77±26 


122.8=b23.2 


20.6=1=1.5 


85.3=bl.5 


5 


HAT-P-34 


-0.004=1= 1.1 64 


4.822±0.179 


24.3±1.2 


0=1=14 


421 ±32 


343.1±21.3 


24.6=1=1.2 


-7=1=12 


6 


WASP- 12 


1.044=1=0.399 


0.640=1=0.210 


1 fi 

1 -°-0.4 


59 

J -20 


204.4=b2.4 


226=b4 


1 d+A/.y 

1 - 4 -0.5 


63 + ^ 


7 


WASP- 16 


0.324=1=0.617 


1.674±0.376 


3.2±0.9 


,,+26 
U -I9 


353±54 


116.7=b2.4 


5.4±0.7 




8 


WASP- 18 


0.774±0.421 


3.234±0.126 


11.2=1=0.6 


13±7 


1768±5 


1818.3=b8.0 


10.9=1=0.6 


12.6±7 


9 


WASP-19 


0.552=1=0.371 


1.984±0.241 


4.4=1=0.9 


15=1=1 1 


200=1=17 


256±5 


3.7=1=0.9 


7=1=13 


10 


WASP-26 


-0.803±0.722 


1.1 17±0.304 


2.2=1=0.7 


■^-26 


137±15 


137.7=1=1.5 


2.6=1=0.9 


-42^ 


11 


WASP-31 


0.076=1=0.130 


2.605=1=0.121 


6.8=1=0.6 


-6±6 


53=1=19 


58.1=1=3.4 


7.2=1=0.6 


-6±3 


12 


G1436 


-0.081±0.357 


-0.192±0.584 


<0.4 




19±3 


18.34±0.52 


<0.4 




13 


Kepler-8 


0.255±0.353 


2.960=1=0.179 


8.9=1=1.0 


5±7 


74±37 


68.4=bl2.0 


8.8=1=0.9 


-5=b6 


14 


HAT-P-2 


0.697±0.341 


4.349±0.101 


19.5=1=1.4 


9±10 


808±9 


983.9=1=17.2 


19.2=1=0.7 


9±4 


15 


HD 149026 


0.558±0.339 


2.706=1=0.157 


7.7=1=0.8 


12±7 


31=1=6 


433=1=1.2 


7.2=1=0.7 


4±6 


16 


HD 209458 


-0.167±0.237 


2.084=1=0.056 


4.4=1=0.2 


-5±7 


128±32 


84.67i0.70 


4.5=1=0.2 


3=b3 


17 



REFERENCES. — (1) Noyes et al. (2008): (2) Pal et al. (2008); (3)Buchhave et al. (2010) ; (4) Kipping et al. (2010); (5) Hartman et al. (2011); (6) Bakos 
et al. (2012); (7) Hebb et al. (2009); (8)Lister et al. (2009); (9) Hellier et al. (2009); (10) Hebb et al. (2010); (1 1) Anderson et al. (201 la); (12) Anderson et al. 
(201 lc); (13) Maness et al. (2007) (14) Jenkins et al. (2010); (15) Pal et al. (2010); (16) Sato et al. (2005); (17) Torres et al. (2008) 

"The results in columns 8-9 are based on an analysis in which a prior constraint was imposed on K+, using the value reported in column 7. We prefer the values 
in columns 4-5 for reasons described in Section 3.5. 

b Reference for the value of K+ in column 7. 

c The internal uncertainty for the A measurement in this system is 14°. The larger uncertainty reported here is the standard deviation of the 3 independent 
measurements, as discussed in section 4.2. 



3.6. Quantitative Analysis 

Having described the model, we now describe the proce- 
dure for parameter estimation. For each observed transit, we 
calculated the expected time of inferior conjunction (T c ) based 
on the published ephemerides. This predicted time was used 
as a prior in our analysis of the transit-night RV data. The out- 
of-transit variation was parameterized by the semi-amplitude 
of the projected stellar reflex motion (K*), and an arbitrary 
additive constant velocity (7). 

The transit geometry was described by the following pa- 
rameters: the radius of the star in units of the orbital semi- 
major axis (R+/a), the radius of the planet in units of the 
stellar radius (R p /R+), and the cosine of the orbital inclina- 
tion (cos/ ). We imposed Gaussian priors on R p /R+, the total 
transit duration from first to last contact (741), and the ingress 
duration (T21), based on previously reported photometric anal- 
yses. By choosing those particular parameters we minimized 
the correlations among their uncertainties (Carter et al. 2008), 
consistent with our treatment of the priors as independent 
Gaussian functions. However, small correlations do exist and 
affect 72i most strongly. In the interest of homogeneity we 
took the simple and conservative approach of doubling the re- 
ported uncertainty in T^i for all systems. Table 2 specifies all 
the priors that were used in practice, including the increased 
uncertainty in T21. 

For systems where values for 741 or T21 were not previously 
reported in the literature we obtained a transit light curve and 
derived those parameters ourselves, using the transit model of 
Mandel & Agol (2002). The uncertainty intervals in R p /R+, 
T41, and 72i were then used as Gaussian priors in the RV anal- 
ysis in the same way as for all the other systems, including the 



doubled uncertainty in T^i. 

We assumed a quadratic limb-darkening law with parame- 
ters mirm and M2RM selected from the tables by Claret (2000) 
for the Johnson V band (similar to the iodine spectral region). 
We used the "jktld" ' tool to query the ATLAS models for all 
systems except Gl 436, for which we used the tables based on 
the Phoenix code. We placed a Gaussian prior on mi RM + «2rm 
with a width of 0.1. The difference wirm - «2rm was held 
fixed at the tabulated value, since this combination is weakly 
constrained by the data and has minimal effect on the other 
parameters. 

Because A and vsim* have correlated uncertainties we used 
the fitting parameters ^/vsini^ cos A and Vvsin/^sin A. A final 
constraint used in our fitting statistic is the prior knowledge 
concerning the projected rotation speed, which was measured 
for all systems in our sample. Here we adopted a minimum 
uncertainty of 1.5 kms" 1 even when the reported uncertainty 
was lower. The vsin/ + measurements in our sample were con- 
ducted by different researchers using different approaches and 
it is not clear that the scales for the different approaches are 
identical, particularly for the very challenging cases of low 
projected rotation speeds. Furthermore we expect that our 
lack of knowledge on differential rotation, macroturbulence, 
and spectroscopic limb darkening introduces systematic er- 
rors at that level. All constraints used in our fits are listed in 
table 2. 

Finally, we must specify the width of a Gaussian function 
representing the contribution to the line width due to both mi- 

1 This code is kindly made available by J. Southworth: 

http : / / www .astro.keele.ac.uk/jkt/codes/jktld. html 
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Fig. 2. — Comparison between results with and without a prior on K*. 

Upper panel: Measurements of Id, based on the transit-night velocity gra- 
dient, versus previously reported measurements of K* from orbital solutions 
based on RV data obtained over months or years. (The three systems with the 
largest K± are omitted to permit the other systems to be viewed more clearly.) 
While for many systems the results are consistent, there are a number of sys- 
tems showing disagreement, possibly due to stellar RV noise with a timescale 
of hours to days. Middle panel: Difference between v sin 4 values obtained 
with and without a constraint on K* based on the orbital solution. Lower 
panel: Same, but for A. 



croturbulence and the PSF of the spectrograph. The influ- 
ence of both can be approximated by a convolution with a sin- 
gle Gaussian, despite their difference in origin (Hirano et al. 
2011c). Assuming a value of 2 kms" 1 for the microturbu- 
lence parameter and 2.2 kms" 1 for the PSF width at 5500 A, 
we obtain a a of 3 km s" 1 for this purpose. The results for 
A and vsini* do not depend strongly on the values of these 
parameters. 

To derive confidence intervals for the parameters we used 
the Markov Chain Monte Carlo algorithm. Before starting 
the chain we also added "stellar jitter" terms in quadrature 
to the internally-estimated uncertainty of the RVs to obtain a 
reduced \~ close to unity. In making this step we assumed that 
the uncertainties in the RV measurements within a given night 
are uncorrelated and Gaussian. Table 5 reports the original, 
internally-estimated uncertainties without any jitter term. The 
added terms are listed for each system in Table 2 

The results reported in Table 3 are the median values of 
the posterior distribution. The quoted uncertainty intervals 
represent the range that excludes 15.85% of the values on each 
side of the posterior distribution, and encompass 68.3% of the 
probability. 

4. NOTES ON INDIVIDUAL SYSTEMS 

In this section we report briefly on each of the 17 systems in 
our sample (14 new observations and three re-analyses of pre- 
viously published data). We mention whenever we deviated 



from the general approach described above. We also discuss 
briefly the values for vsin/* and A that were obtained. For 
cases in which A has been previously reported in the litera- 
ture, we compare the new value with the previous value. 

4.1. HAT-P-6 

The discovery of the hot Jupiter HAT-P-6 b was reported by 
Noyes et al. (2008). We obtained 42 HIRES spectra over a 
five-hour period covering the transit, which occurred on the 
night of 5/6 November 2010. Although little data could be 
obtained before ingress, the observations continued for 1.5 
hours after the 3.5-hour transit. All the RVs are shown in the 
left panels of Figure 3 along with the best-fitting model. The 
RV curve immediately reveals a retrograde orbit. The right 
panels show the results in the A - vsin/ + plane. To obtain these 
we used the HIRES RVs together with the prior information 
from Noyes et al. (2008) which are reported in Table 2. We 
find A= 165±6° and vsim* = 7.8 ± 0.6 kms" 1 . With this 
we confirm, with higher precision, the measurement presented 
by Hebrard et al. (2011a). They found A = 166 ± 10° and 
vsin/* = 7.5 ± 1.6 kms" 1 . 

When the model was rerun with a prior constraint on X*, 
we obtained A = 175 ±4°, in disagreement with the prior-free 
result and with the previously reported value. (See Table 3 
and Figure 2.) Apparently the star exhibits RV noise that is 
correlated on the timescale of at least a few hours. 

4.2. HAT-P-7 

Pal et al. (2008) reported the discovery of HAT-P-7b. This 
was one of the first systems for which a significant misalign- 
ment between the stellar spin and orbital angular momentum 
was discovered in a star-planet system (Winn et al. 2009; 
Narita et al. 2009a). Both research groups reported an appar- 
ently retrograde orbit. In fact this star is likely to be viewed 
nearly pole-on, judging from its unusually low vsin/* (Winn 
et al. 2009; Schlaufman 2010), but with nearly antiparallel sky 
projections of the orbital and stellar rotation vectors. 

Interestingly the two independent results for A were in dis- 
agreement, with Winn et al. (2009) reporting A of 182.5 ±9° 
and vsin/ + = 4.9^1 kms" 1 , and Narita et al. (2009a) re- 
porting A = -132.6+}^° (equivalent to A = 227.4^j°) and 
vsin/* = 2.3^} 5 km s" 1 . Both groups used the same instrument 
(HDS on Subaru). With HIRES we reobserved the system be- 
fore, during, and after a planetary transit occurring during the 
night 23/24 July 2010. The RVs are displayed in Figure 4. 
While the RM effect is clearly detected, the SNR of the de- 
tection is low enough that an analysis of the RM effect would 
benefit from a very precise prediction for the transit mid- 
point. For this reason we downloaded the Quarter 6 Kepler 
light curve for this system, covering the epoch of our spec- 
troscopic transit observations. The analysis of transit light 
curves is described in section 3.6 and the derived photometric 
priors are given in Table 2. We obtain a A of 155 ± 14° and a 
vsin/* = 2.7 ±0.5 kms" 1 . 

While the v sin /* is intermediate between the two values re- 
ported previously, and the finding of a retrograde orbit is con- 
firmed, the value we obtained for A does not agree with either 
of the previously reported values. What can have caused the 
disagreement between the three results? While we can not 
give a definite answer we suspect that two characteristics of 
this system might have been important. K+ is ten times larger 
than the amplitude of the RM effect making a clean separation 
of the RM effect from the orbital radial velocity challenging. 
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FIG. 3 . — Spectroscopy of HAT-P-6 transit Left panels- The radial velocities measured before, during, and after transit are plotted as a function of time from 
inferior conjunction. The black error bars indicate the internal RV uncertainties. The gray error bars also include "stellar jitter" as explained in the text. The upper 
panel shows the measured RVs, with 7 subtracted, and the best-fitting model. The lower panel shows the RVs after subtracting the best-fitting orbital model. 
The light and dark gray bars in the lower panel indicate times of first, second, third, and fourth contact. Right panels- The gray scale indicates the posterior 
probability density, marginalized over all other parameters. The black, dark gray, and light gray contours represent the 2-D 68.3%, 95%, and 99.73% confidence 
limits. The one-dimensional marginalized distributions are shown on the sides of the contour plot. The dash-dotted line shows the prior applied to vsini*. 
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FIG. 4. — Spectroscopy of HAT-P-7 transit Similar to Figure 3 but this time for data obtained during a transit in the HAT-P-7 system. In the middle panel on 
the left side, the orbital contribution to the observed RVs has been subtracted, isolating the RM effect. 



In addition the RM amplitude is with ~ 20 ms" 1 only a few 
times larger than the typical uncertainty in the RV measure- 
ments. Albrecht et al. (2011a) argued that for low SNR de- 
tections and for A near 0° or 1 80° the uncertainty in A can be 
underestimated. However we cannot rule out the possibility 
that our model is missing some aspect of stellar astrophysics, 
perhaps one that is most apparent for nearly pole-on systems. 
The HAT-P-7 system is the only system in our sample where 
there exists a significant disagreement between several inde- 
pendent measures of A. 



Since the discrepancies are unresolved, and to avoid over- 
interpretation of the results in the subsequent discussion, the 
uncertainty reported in Table 3 is the standard deviation (37°) 
of the three independent measurements. 

4.3. HAT-P-16 

The existence of HAT-P-16b was announced by Buchhave 
et al. (2010). We observed this northern system with the 
Keck I telescope during the night 26/27 December 2010. Due 
to bad weather we were only able to observe the second half 
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FIG. 6. — Spectroscopy of HAT-P-24 transit Similar to Figure 3 but this time for data obtained during a transit in the HAT-P-24 system. 



of the transit (Figure 5). No data before the transit were ob- 
tained. Using these RVs together with the photometric and 
vsin/ + (3.5 ±1.5 kms" 1 ) priors from Buchhave et al. (2010) 
we found that HAT-P-16b orbits its host star on a prograde 
orbit (A = 2+46°). While the data are consistent with good 
alignment, little more can be learned. Concerning vsim* we 
found a value of 3.1 ± 1.0 kms" 1 , consistent with the prior. 
Moutou et al. (2011) found a projected obliquity of 10 ± 16° 
for this system. In the subsequent discussion of Section 5, we 
adopt their value. 

4.4. HAT-P-24 

The 24th HATNet planet was announced by Kipping et al. 
(2010). We obtained HIRES spectra of this system during the 



night 27/28 September 2010. As was the case for HAT-P-16, 
the observations were interrupted by bad weather, but at least 
for HAT-P-24 we obtained one pre-ingress data point (see Fig- 
ure 6). This is why the result, A = 20 ± 16°, is much more 
precise than for HAT-P-16. Nevertheless with the current data 
set is it not possible to exclude a small misalignment. 

4.5. HAT-P-32 

HAT-P-32 b was detected around a star which displays RV 
jitter on the order of 80 m s" 1 (Hartman et al. 201 1). Therefore 
the SNR of the RM effect detection in our dataset, obtained 
on 5/6 December 201 1, is relatively low (see Figure 7). Nev- 
ertheless this is the system in our sample for which we ob- 
tain the highest precision in our measurement of the projected 
obliquity (A = 85 ± 1 .5°). This curious fact is the consequence 
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FIG. 7. — Spectroscopy of HAT-P-32 transit Similar to Figure 3 but this time for data obtained during a transit in the HAT-P-32 system. 
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Spectroscopy of HAT-P-34 transit Similar to Figure 3 but this time for data obtained during a transit in the HAT-P-34 system. 



of the system having a low impact parameter (0. 12) and at the 
same time presenting an asymmetric RM curve. 

For systems with low impact parameters, even a strong mis- 
alignment will lead mostly to a reduction of the amplitude of 
the RM effect, with hardly any asymmetry in the RM curve, 
(e.g. Gaudi & Winn 2007; Albrecht et al. 201 lb). Only a nar- 
row range of angles near 90° can produce a asymmetric RM 
effect. Therefore when an asymmetry is detected in such a 
system, even relatively coarse RV data will lead to high pre- 
cision in A. Interestingly if we take the vsin/* prior at face 
value and assume that our RM model contains all the relevant 
physics, then we find with T 2 \ = 0.41251 ± 0.0017 hr a five- 
fold improvement in our knowledge of the ingress duration 
compared to the prior constraint (72.] = 0.4128 ± 0.0096 hr). 



This is because for a given vsim* and A near ±90°, the im- 
pact parameter is encoded in the amplitude of the RM effect. 

The pairing of HAT-P-32 and HAT-P-2 is instructive (see 
Section 4. 15). For HAT-P-2 the RM effect was observed with 
a very high SNR: the ratio between the amplitude of the RM 
effect and the typical RV uncertainty in the RVs is ~ 10. But 
despite the very high SNR in the detection of the RM effect, 
we obtained a relatively low precision in the measurement of 
A. This is the reverse of what we see for HAT-P-32. For this 
system, the ratio between the amplitude of the RM effect and 
the average RV uncertainty is ~ 2. Nevertheless we obtain a 
very high precision in the measurement of A. This highlights 
how important the geometry of the transit together with the 
projected obliquity is in determining the pression in the mea- 
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FIG. 10. — Spectroscopy of WASP-16 transit Similar to Figure 4 but this time for data obtained during a transit in the WASP 16 system. 



surement of the projected obliquity. 

4.6. HAT-P-34 

HAT-P-34b was announced by Bakos et al. (2012) together 
with three other transiting planets. HAT-P-34 b orbits its host 
star on a highly eccentric (<? = 0.44) orbit with a relatively long 
period of 5.5 days. We observed this system with HIRES dur- 
ing the night 2/3 September 2011. We obtained 62 spectra 
during and after the transit (Figure 8). As the RM curve sug- 
gests we find that the projections of the stellar rotation and 
orbital angular momenta are well aligned (A = 0± 14°). 

4.7. WASP-12 

WASP- 12 harbors a transiting exoplanet with a short orbital 
period of It was discovered by Hebb et al. (2009). We 



took from those authors our prior on vsin/*, while we ob- 
tained from Maciejewski et al. (2011) our priors on the geo- 
metric transit parameters. 2 

We observed the RM effect in this system on 1/2 Jan- 
uary 2012 with the Keck I telescope. Our analysis of the 38 
RVs (Figure 9, left panel) obtained before, during, and af- 
ter the 3 hr transit indicates a misalignment of 59^0° an d 
vsin/* = 1.6^0 4 kms -1 . For this system there is a strong cor- 
relation between A and vsin/*; for higher projected rotation 
speeds A increases and approaches 90° (Figure 9, right panel). 
This is interesting as WASP-12, with a mass of 1 .35 M Q (Ma- 
ciejewski et al. 2011), is expected to be a fast rotator and to 

2 Husnoo et al. (201 1) also obtained RV data during a transit of this system, 
but their results were not conclusive. 
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have a high v. Using the method employed by Schlaufman 
(2010) we find an expected v of 13.7 ±2.5 kms" 1 suggesting 
a low sin/*, i.e., the stellar spin axis is inclined along the line 
of sight. It seems very likely that the star and the planetary or- 
bit are misaligned, though we cannot tell with certainty how 
much of the misalignment is in the plane of the sky as op- 
posed to the perpendicular direction. A similar situation was 
found for WASP-1 (Albrecht et al. 201 lb). 

4.8. WASP-16 

Lister et al. (2009) found this transiting hot Jupiter on a 
3.1 -day orbit around a southern-sky star. We used the PFS 
in conjunction with the Magellan II telescope to obtain RVs 
during a transit occurring 3/4 June 2010 (Figure 10). We used 
the ephemeris and information on the projected rotation speed 



from the discovery paper as prior information in the fit. As no 
information on T2\ was given, we used the EULERCAM light 
curve presented in the discovery paper to establish priors for 
Til, T21, and R p /R+. This system is one of those for which 
we find (based on the fit to the out-of-transit observations) an 
orbital velocity semi-amplitude (K\ = 353 ±54 ms" 1 ) that is 
significantly different from the orbital solution presented in 
the discovery paper (Ki = 1 16.7+j g m s -1 ). To investigate this 
issue we obtained 10 out-of-transit observations on a number 
of different nights. These observations agree with the previ- 
ously reported orbital solution. A large stellar spot might be 
responsible for the excess RV change during the transit night 
(see section 3.5), although such a spot would need to cover a 
substantial part of the stellar surface. 
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FIG. 13. — Spectroscopy of WASP-26 transit Similar to Figure 4 but this time for data obtained during a transit in the WASP 26 system. 
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FIG. 14. — Spectroscopy of WASP-31 transit Similar to Figure 3 but this time for data obtained during a transit in the WASP31 system. 



We omitted 9 data points obtained at the end of the tran- 
sit night, after the template observation for WASP- 16 was 
obtained. These observations would in principle facilitate a 
separation of the RM signal from other sources of RV vari- 
ation. However, given that there is already an indication of 
intrinsic RV variations apart from orbital motion, and the ab- 
sence of any guarantee that the intrinsic variations are lin- 
ear in time over the course of many hours, we decided to 
use only the data points immediately following egress. All 
the RVs are listed in Table 5. Using the 31 remaining RV 
data points obtained during and after the transit we obtained 
A= 11^° and v sin/, = 3. 2 ± 0.9 kms" 1 . Brown et al. (2012) 
recently measured A = -4. 2^3 9 ° consistent with the value ob- 
tained here. However they found a significantly lower value 



for vsin;* (1.2^'j kms" 1 ) which is not consistent with their 
spectroscopic prior, nor with ours. We also note that Brown 
et al. (2012) did not report any disagreement between the out- 
of-transit velocity gradient observed on the transit night, and 
the published orbital solution. 

4.9. WASP- 18 

WASP-18b orbits its host star in only 0.94 days (Hellier 
et al. 2009). As noted by the discoverers it is an excellent 
system to study the tidal interaction between a gas giant and 
its host star. Triaud et al. (2010) found the projections of the 
stellar spin and orbital axes to be aligned (A = 4 ±5°). 

We obtained 48 data points with the PFS during the night 
8/9 October 201 1. Six of these data points were obtained at a 
significant distance in orbital phase from the transit, and are 
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FIG . 15. — Spectroscopy of Gl 436 transit Similar to Figure 3 but this time for data obtained during a transit in the GI 436 system. 
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FIG. 16. — Spectroscopy of Kepler-8 transit Similar to Figure 3 but this time for data obtained during a transit in the Kepler 8 system. Black filled circles 
represent the RV data points already obtained by Jenkins et al. (2010), while the open circles show the new HIRES data. 




not included in the fit, for reasons discussed in the preceding 
section and section 3.5. Figure 1 1 presents the remaining 42 
RVs obtained before, during, and after the transit. 

No value for T^i has been reported in the literature for 
WASP- 18 b. We obtained the photometry from South worth 
et al. (2009) and performed our own analysis. For the rea- 
sons mentioned by Southworth et al. (2010) we did not use 
this light curve to improve on the ephemeris for this system. 
The photometric priors reported in Table 2 for this system are 
derived from the fit to that light curve. 

We found that the projected stellar spin (A = 13 ±7°) ap- 
pears to be misaligned with the orbital rotation axis, although 
the statistical significance is marginal. The value for vsin/* 
11.0 ±0.5 kms" 1 does not give any strong indication of an 



inclination along the line of sight. Any obliquity in WASP-18 
seems to be small. 

4.10. WASP-19 

WASP- 19 is another system with a very short-period 
(0.78 d) hot Jupiter discovered by the WASP consortium 
(Hebb et al. 2010). We observed this system with PFS during 
the transit night 20/21 May 2010 (Figure 12). Our photomet- 
ric priors were based on results by Hellier et al. (201 1). These 
authors derived the projected obliquity in the same fit that was 
used to determine the photometric parameters. This could in 
principle make our observation of A dependent on their obliq- 
uity observation. However given the quality of their photo- 
metric and RV data, we expect that the RVs obtained by them 
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FIG. 17. — Spectroscopy of HAT-P-2 transit Similar to Figure 4 but this time for data obtained during a transit in the HAT-P-2 system. Structure can be seen 
in the residuals indicating that our modeling does not capture all the relevant physics of the RM effect. 
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FIG. 18. — Spectroscopy of HD 149026 transit Similar to Figure A but this time for data obtained during a transit in the HD 149026 system. 



have little influence on their results for T^\, Ti\, and R p /R+. 
We therefore expect our result (A = 15 ± 1 1°) to be indepen- 
dent of their result (A = 4.6 ± 5.2°). Our prior on the projected 
rotation speed was taken from the discovery paper. 

As can be seen in Figure 1 2 there is correlated excess noise 
in the RVs during pre-egress. This together with the low SNR 
detection raises the concern that we might not have detected 
a RM effect at all, and we are fitting our RM model to noise. 
We carried out an experiment similar to the one described by 
Albrecht et al. (2011b) for the WASP-2 system. Briefly, we 
create 200000 fake datasets with similar noise characteristics 
as the real dataset, but without an underlying RM effect. We 
then fitted an RM model to each dataset. Form the result- 
ing density distribution we conclude that our RM detection 



represents a 3a detection. We also used the Bayesian Infor- 
mation Criterion (BIC) to assess whether the model including 
the RM effect is preferred over a model without the RM effect 
(e.g. Brown et al. 2012). The RM model is preferred, with 
BICrm = 48 and BIC no 

rm — 60. However the BIC test de- 
pends somewhat on how many of the RV measurements were 
obtained during transit, compared to the number obtained out- 
side of transit. Furthermore both tests (or at least our im- 
plementations of these tests) assume that the noise is uncor- 
rected, which does not seem to be the case for the present 
dataset. 



4.11. WASP -26 
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FIG. 19. — Spectroscopy of HD 209458 transit Similar to Figure 4 but this time for data obtained during a transit in the HD 209458 system. 



This system was detected by Smalley et al. (2010). An- 
derson et al. (2011a) attempted a measurement of the RM 
effect but did not make a secure detection. We obtained 39 
RV data points during the night 17/18 August 2010 (Fig- 
ure 13). Using these RVs, photometric priors from Ander- 
son et al. (201 la), and a prior on the projected stellar rotation 
speed from Smalley et al. (2010) 1 we derive A = —34^6° and 
vsin/i = 2.2 ±0.7 kms" 1 . 

4.12. WASPS 1 

The existence of WASP-3 1 b was announced by Anderson 
et al. (2011c). We observed this system with the Keck I tele- 
scope during the night 12/13 March 2012 and obtained 27 RV 
measurements (Figure 14). 

Using these RVs together with the vsin/*, Z41, T%\, and 
Rp/R* priors from Anderson et al. (2011c) and priors on 
the ephemeris from Dragomir et al. (2011) we found a low 
projected obliquity (A = -6 ± 3°). We obtained vsim* = 
7.3 ±0.4 kms -1 , consistent with the prior (7.9 ±1.5 kms" 1 ). 
Brown et al. (2012) recently found A = 2.8±3.1° and v sin z* = 
7.5 ±0.7 kms -1 . The value for A is inconsistent with our 
result. Brown et al. (2012) used the ephemeris presented in 
Anderson et al. (201 lc). Using the T c and P values from An- 
derson et al. (2011c) to set the constraint on the time of in- 
ferior conjunction for the observed spectroscopic transit we 
obtain A = 2 ± 3°. This is consistent with the result by Brown 
et al. (2012). For our final results we decided to use the tim- 
ing information from Dragomir et al. (2011), who used the 
ephemeris from the discovery paper in combination with an 
additional light curve (although we note that the new light 
curve contains only a small amount of post-egress data). 

Therefore while we have excellent agreement between two 
different research groups, the dependence of the result on 
the photometric priors and the inconsistency between the re- 

3 We used the vsini* value from Smalley et al. (2010) as prior, and not the 
value reported by Anderson et al. (201 la). This is because the former authors 
estimated v sin i„ under the assumption that the macroturbulent parameters is 
4.1 ± 0.3 km s . This is closer to our estimation (5.1 ± 1.5 km s ) than the 
value of 3.0 ±0.3 km s" 1 used by the latter authors. 



suits using the same RVs but different timing information cast 
some doubt on the formal uncertainty intervals. For this rea- 
son, we double the uncertainty in A before including this sys- 
tem in the subsequent discussion of the interpretation of all 
the results. Future photometric observations will be helpful, 
but at least it seems clear that the projected obliquity is small. 

4.13. G1436 

The transiting planet in G1436 was discovered by Butler 
et al. (2004) and found to be transiting by Gillon et al. (2007). 
It would be of interest to know the obliquity in this system as 
G1436b is of similar mass then HAT-P-llb, which is so far 
the only Neptune-class planet for which the host star's obliq- 
uity has been measured. Using HIRES we obtained RVs of 
the system during the night of 24/25 April 2010. To reduce the 
uncertainties in the photometric parameters, we gathered new 
photometric data with Keplercam, a CCD camera on the 1.2 m 
telescope of the Fred L. Whipple Observatory on Mount Hop- 
kins, Arizona (Szentgyorgyi et al. 2005). The RV measure- 
ments are displayed in the left panels of Figure 15. The right 
panels display the posterior in the A - vsin/* plane. No RM 
effect was detected, and therefore the data provide no con- 
straint on A. Only an upper limit of 0.4 kms" 1 on vsim* can 
be obtained. This upper limit on vsin/* is more constraining 
than the prior we adopted, which was a one-sided Gaussian 
prior enforcing vsinz* < 5 kms" 1 . Thus from the transit data 
we learned only that the star is a slow rotator. The structure in 
the posterior distribution is caused by the lower sensitivity of 
our measurement for nearly prograde and retrograde configu- 
rations. See Albrecht et al. (201 lb) for an explanation of this 
effect. 

4.14. Kepler-8 

Jenkins et al. (2010) reported the discovery of Kepler-8 and 
also found A = -26 ± 10.1°, a misalignment between the pro- 
jections of the stellar and orbital spins. As the data set avail- 
able to Jenkins et al. (2010) only consisted of data taken dur- 
ing the transit itself, it is difficult to assign a significance to the 
result. We therefore obtained a new data set with HIRES dur- 
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ing the transit night 7/8 August 2011. We used both HIRES 
data sets in our analysis, including two additional parame- 
ters representing a second velocity offset and a second RV 
transit-night velocity gradient between the two different tran- 
sit nights. Flat priors were assigned to these parameters. We 
found that the data are consistent with good alignment be- 
tween the stellar and orbital rotation axes (see Figure 16). 
We obtained A = 5 ± 7° and vsim* = 8.9 ± 1.0 kms" 1 . The 
velocity offset between the two data sets was found to be 
12 ± 10 ms" 1 . Using only the new dataset in our analysis, 
we obtained A = 5 ± 10°. Using only the older dataset, we 
obtained A = 30^|i °- We chose to use both datasets for this 
system, as they have been obtained with the same instrument 
and setup, and the RVs were extracted with the same code, 
circumstances different from all the other systems for which 
multiple transits have been observed. Furthermore this case 
highlights the importance of prior constraints on the data. 
Without a constraint on the RV offset as applied by Jenkins 
et al. (2010), the original dataset is consistent with good align- 
ment. 

4.15. HAT-P-2 

HAT-P-2 was one of the first systems for which a measure- 
ment of the projected obliquity was reported. Winn et al. 
(2007) measured A = 1.2 ± 13.4°. Loeillet et al. (2008) found 
A = 0.2^25°, a similar value. The uncertainty in this measure- 
ment is relatively high because HAT-P-2 b has a low impact 
parameter. After these results were published, new and im- 
proved system parameters based on new photometry were re- 
ported by Pal et al. (2010). For this reason we reanalyzed the 
RVs obtained by Winn et al. (2007) together with the new pho- 
tometric priors, obtaining A = 9 ± 5°. Thanks to the improved 
photometric information, the formal uncertainty in A is now 
lower, and supports a 1.8cr "detection" of a misalignment, a 
different conclusion from the one by Winn et al. (2007). How- 
ever, one can see in the lower left panel of Figure 1 7 that, after 
subtraction of our best fitting model, structure remains in the 
residuals. 

This is a clear sign that our model does not capture all the 
effects influencing the RV anomaly during transit. We there- 
fore consider the uncertainty of our analysis to be itself quite 
uncertain, by at least a factor of two. In Table 3 we report 
a doubled uncertainty relative to our formal results for A and 
vsin/ + . Therefore while the photometric information as well 
as the quality of the RV data would allow for a measurement 
of A to within a few degrees, the limitations of our model 
of the RM effect prevent this gain from being fully realized. 
The model expects a larger RM effect shortly before and after 
mid-transit and a lower amplitude around the second and third 
contacts. One possible effect that would produce this type of 
structured residuals is a change of limb darkening with depths 
of stellar absorption lines. We note that no such structure in 
the residuals was seen in the original analysis of these data 
using a different modeling for the RM effect, which was de- 
rived empirically for the specific spectrograph and stellar type 
(Winn et al. 2007). 

It is not surprising that the limitations of our model would 
be most apparent for HAT-P-2, as it has the combination of 
very rapid rotation and very high signal-to-noise ratio, given 
the brightness of the star (V = 8.7). For this system it would 
be better to analyze the absorption line profiles directly, and 
their changes over the course of the night, rather than deriving 
and modeling RVs (e.g. Albrecht et al. 2007; Collier Cameron 



et al. 2010). 

4.16. HD 149026 

Wolf et al. (2007) found A = -12 ± 15° for this well-known 
system. Since that time, (Carter et al. 2009) presented a more 
precise light curve based on Hubble Space Telescope (HST) 
infrared observations. Thus a re-analysis is warranted. We 
used the photometric information from Carter et al. (2009) 
and the vsim* prior from the discovery paper (Sato et al. 
2005). We found a projected obliquity of 12 ±7°. This is 
only marginally consistent with the previously reported value 
of A using the same RV data. The main reasons for the dif- 
ference are the use of new photometric data, and the lack of 
post-egress data (Figure 18) in combination with our choice 
to not impose a prior constraint on K*. When we repeated the 
analysis with a prior on (Table 3) we found the results to 
be consistent with those of Wolf et al. (2007). 

4.17. HD 209458 

The first system in which a transiting planet was discovered 
was also the first system for which the projected obliquity was 
measured. Queloz et al. (2000) found an angle consistent with 
alignment to within 20°. Winn et al. (2005) found A to be 
-4.4 ±1.4°. The later result suggested that a small but signif- 
icant misalignment exists in this system, similar to the mis- 
alignment of our Sun against the ecliptic (~ 7°). However, 
Winn et al. (2005) used data from different nights, unlike our 
current procedure. If the star is active or has star spots, the 
long-term intrinsic RV noise would be manifested as velocity 
offsets between different nights, which were not taken into ac- 
count by Winn et al. (2005). Indeed we find evidence that the 
stellar jitter has different amplitudes on different timescales: 
by fitting a model to all of the available HIRES data, we find 
the rms residual of the data obtained on different nights to be 
5 m s -1 while the rms residual of the transit-night data is only 
3 ms" 1 . 

We therefore repeated the analysis using only the HIRES 
RVs obtained during the transit night, 29/30 July 2000. We 
also used the HST light curve obtained by Brown et al. (2001) 
to obtain photometric priors, and used the vsin/* value from 
Laughlin et al. (2005) as a prior. We obtained a projected 
stellar obliquity of -5 ±7°. The looser bounds on A are a 
consequence of the lack of egress and post-egress data on the 
night 29/30 July 2000 (Figure 19). This case emphasizes the 
need to obtain data outside of the transit, unless one is willing 
to assume that the long-term RV noise is negligible. 

5. DISCUSSION 

We now use the measurements of projected obliquities ob- 
tained in the last section together with other measurements 
taken from the literature to learn more about the evolution of 
these systems. 4 

In 5. 1 we first remark on some specific cases from the liter- 
ature. In Section 5.2 we revisit the evidence for an correlation 
between the degree of alignment and the effective stellar tem- 
perature discovered by Winn et al. (2010a). We then discuss 
further evidence of tidal interaction, based on correlations be- 
tween the planet-to-star mass ratio and A, and the dependence 
of A on the orbital distance between the two bodies. In Sec- 
tion 5.2 we sort the systems according to a theoretical tidal 

4 We present here all measurements with the sign convection form Ohta 
et al. (2005) and use the symbol A. Often researchers use the symbol (5 and 
the sign convention from Hosokawa (1953), where A = — j3. 
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TABLE 4 
Measured projected obliquities 





System 


(km s _1 ) 


\ 

C°1 
i. ) 


ret 


1 


CoRoT-2 b 


10.9±0.5 


4-OjJ;? 


1 


2 
3 


CoRoT-3 b 
CoRoT-18b 


17.0±1.0 
8.0±1.0 


37.6^3 
-10±20 


2,3 
4 


4 


HAT-P-1 b 


3.8±0.6 


3.7±2.1 


5 


5 


HAT-P-2 b 


19.5±1.4 


9±10 


this work 


6 


HAT-P-4 b 


5.8±0.3 


-4.9±11.9 


6 


7 


HAT-P-6 b 


7.8±0.6 


165±6 


this work 


8 


HAT-P-7 b 


2.7±0.5 


155±37 


this work 


9 


HAT-P-8 b 


14.5±5 


-17+ 92 

-11 5 


7 


10 


HAT-P-9 b 


12.5±1.8 


-16±8 


7 


11 


HAT-P-1 1 b 


1.0±0.9 


103^6 


8 


12 


HAT-P-13b 


1.7±0.4 


1.9±8.6 


9 


13 


HAT-P-14b 


8.2±0.5 


-170.9±5.1 


6 


14 


HAT-P-16b 


3.9±0.8 


-10±16 


7 


15 


HAT-P-23 b 


7.8±1.6 


15±22 


7 


16 


HAT-P-24 b 


11.2±0.9 


20±16 


this work 


17 


HAT-P-30 b 


3.1±0.2 


73.5±9.0 


10 


18 


HAT-P-32 b 


20.6±1.5 


85±1.5 


this work 


19 


HAT-P-34 b 


24.3±1.2 


0±14 


this work 


20 


HD 17156 b 


4.1±0.3 


10±5.1 


11 


21 


HD 80606 b 


1.7±0.3 


42±8 


12 


22 


HD 149026 b 


7.7±0.8 


12±7 


this work 


23 


HD 189733 b 


3.1±0.1 


-0.5±0.4 


13 


24 


HD 209458 b 


4.4±0.2 


-5±7 


this work 


25 


KOI- 13 b 


65±10 


24/154±4 


14 


26 


Kepler-8 b 


8.9±1.0 


5±7 


this work 


27 


Kepler- 17 b 


4.7±1.0 


< 15 


15 


28 


TrES-1 b 


1.1±0.3 


30±21 


16 


29 


TrES-2 b 


2±1 


-9±12 


17 


30 


TrES-4 b 


8.5±1.2 


6.3±4.7 


18 


31 


WASP-1 b 


0.7±1.4 


-59±99 


-26 


32 


WASP-3 b 


14.1±1.5 


3 3 +2 ' 5 


20 


33 


WASP-4 b 


2.1±0.4 


i+14 
1 -12 


21,22 


34 


WASP-5 b 


3.2±0.3 


12 1+ 8 


22 


35 


WASP-6 b 


1.6±0.3 


_11+18 
1J -14 


23 


36 


WASP-7 b 


14±2 


86±8 


24 


37 
38 
39 


WASP-8 b 
WASP- 12 b 
WASP- 14 b 


20±0.6 

1 6 +os 

1 - -0.4 

2.8±0.6 


-123+ 3 - 4 

J -20 

-33.1±7.4 


25 

this work 
26 


40 


WASP- 15 b 


4.3±0.4 


-139.6+4| 


22 


41 


WASP- 16 b 


3.2±0.9 


ir 26 ' 

-19 


this work 


42 
43 


WASP- 17 b 
WASP-18b 


9.9±0.5 
11.2±0.6 


-148.7^7 
13±7 


27 

this work 


44 


WASP- 19 b 


4.4±0.9 


15±11 


this work 


45 


WASP-22 b 


4.4±0.3 


22±16 


28 


46 


WASP-24 b 


7±0.9 


-4.7±4 


29 


47 


WASP-25 b 


7±0.9 


14.6±6.7 


30 


48 


WASP-26 b 


2.2±0.7 




this work 


49 


WASP-3 1 b 


6.8±0.6 


-6±6 


this work 


50 


WASP-33 b 


86±1 


-107.7±1.6 


31 


51 


WASP-38 b 


8.6±0.4 


15+ 33 


29 


52 


XO-3b 


17.0±1.2 


37.3±3.0 


32 


53 


XO-4b 


8.8±0.5 


-46.7±8.1 


33 
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timescale, as a test of whether tides have been important in 
altering the obliquities of these systems. We also discuss the 
implications for the origin of hot Jupiters and the strengths 
and weaknesses of our interpretation. 

5.1. Remarks on the previous literature 

CoRoT-1 — The RM effect was measured by Pont et al. 
(2010). We do not include this measurement in our subse- 
quent discussion because of the low SNR of the detection and 
because it was only possible to obtain a few out-of-transit ob- 
servations. The authors caution that systematic uncertainties 
could cause the actual errors in the measurement of the pro- 
jected obliquity (A = 77 ± 1 1°) to be larger than the statistical 
uncertainty. We note that if the value for A from this study 
were taken at face value, it would constitute an exception to 
the pattern presented below. For this reason it would be an im- 
portant system to reobserve. It is a challenging target because 
of the faint and slowly-rotating host star. 

CoRoT-11 — Gandolfi et al. (2010) obtained RVs during a 
planetary transit in this system. As only the first half of the 
transit was observed, and only with a low SNR, they could 
only conclude that the orbit was prograde. In that sense this 
is a similar case to HAT-P-1 6 presented here. We decided to 
exclude CoRoT-1 1 from subsequent discussion as we did with 
our result on HAT-P- 16. 

CoRoT-19 — This system has an F9V star and a Jupiter-mass 
planet on a 3.9-day orbit. The RM effect was measured by 
Guenther et al. (2012). They found A = -52^2 °> a misalign- 
ment between stellar and orbital axes. However the RM effect 
was only detected at the 2. 3c level, and no post-egress data 
were obtained. We omit this measurement in the subsequent 
discussion. 

KOI-13 — Szabo et al. (2011) detected a slight asymmetry 
in the transit light curve and attributed the asymmetry to a 
misalignment of the planet's orbit relative to the stellar rota- 
tion. The host star is a fast rotator, leading to a lower surface 
gravity and surface brightness around the equator compared 
to the poles. By modeling this effect, Barnes et al. (2011) 
found |A| to be either 24 ±4° or 156 ±4°. Either choice rep- 
resents a substantial misalignment and would lead to similar 
conclusions in the subsequent discussion. For simplicity of 
presentation in the plots to follow, we arbitrarily adopt the 
lower value. Barnes et al. (2011) also calculated the stel- 
lar inclination along the line of sight, which we do not use 
here, since this information is not available for most of the 
other systems. We adopt the mass for the secondary estimated 
by Mislis & Hodgkin (2012), which is also consistent with 
work by Shporer et al. (2011). Both estimates were based on 
the photometric orbit. We further adopt the age estimate of 
(710!|^ M y r ) b Y Szab ° et aL (2011). 

WASP-23 — We omit the measurement by (Triaud et al. 
2011). Because of a low impact parameter, the only conclu- 
sion that could be drawn is that orbit is prograde. 

XO-2 — We omit the measurement by (Narita et al. 2011). 
They found A = 10 ±72°, a prograde orbit but with a very 
large uncertainty, similar to our result for HAT-P-16. 

5.2. Relevant system properties 

Effective temperature — Winn et al. (2010a) noted that for hot 
Jupiters, nonzero values of A tended to be associated with hot 
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Fig. 20. — Projected obliquities and projected stellar rotation speeds as a function of the stellar effective temperature. Upper panel: Measurements 
of projected obliquities as a function of the effective temperature of the host star. Stars which have temperatures higher then 6250 K are shown with red filled 
symbols. Blue open symbols show stars with temperatures lower then 6250 K. Stars which measured effective temperature include 6250 K in their 1-cr interval 
are shown by split symbols. Systems which harbor planets with mass < 0.2 Mj up or have an orbital period more then 7 days are marked by a black filled circle 
with a ring. Lower panel: Projected stellar rotation speeds v sin u of the stars in our sample. In addition v sin i„ measurements of stars in the catalogue by Valenti 
& Fischer (2005) are shown as small dots. 



stars. For effective temperatures > 6250 K the obliquities 
have a broad distribution, while for lower temperatures the 
measurements are consistent with low obliquities. The only 
exceptions were two systems with significantly longer orbital 
periods than the rest. Schlaufman (2010) independently found 
that hot stars tended to be misaligned with the orbits of hot 
Jupiters by comparing the measured value of vsku* with the 
expected value for v, for a star of the given mass and age. 
While this approach has the virtue of requiring less intensive 
observations, it does rely heavily on accurate measurement 
of vsin/^, which is problematic for slowly -rotating stars, and 
assumes that v was not affected by any tidal influence of the 
close-in planet. 

With the new measurements presented in this paper and 
with measurements by others over the last two years, the num- 
ber of systems with measured projected obliquities is now up 
to 53. This is more than twice the number that was available 
to Winn et al. (2010a). Table 4 shows the measured projected 



obliquities of all systems used in this study. 1 The increased 
number in measurements enables a stringent test of the pro- 
posed pattern, as well as a more in depth analysis and com- 
parison to other system parameters. 

It should be noted that in almost all cases, the measure- 
ments are of the projected obliquity, and not the true obliquity. 
For true obliquities smaller then 90° the projected obliquity is 
usually smaller than the true obliquity, while for true obliq- 
uities >90° the projected obliquity is usually larger than the 
true obliquity. This factor complicates any detailed compar- 
ison between the measurements and the theoretical expecta- 
tions. For simplicity we have chosen to work with projected 
obliquities, rather than attempting any deprojection scheme 
(Fabrycky & Winn 2009; Morton & Johnson 2011). 

The upper panel of Figure 20 shows the projected obliqui- 
ties plotted as a function of the effective temperature of the 

5 See also exoplanets.org, exoplanet . eu, www. astro, 
keele . ac . uk/ jkt/tepcat, and www.aip.de/People/rheller 
for listing of obliquity measurements. 
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companion the degree of misalignment decreases, while for some of the more massive planets a small but significant obliquity is detected. These planets orbit 
stars with radiative envelopes, reducing the effectiveness with which tidal energy can be dissipated. Names of systems with particular small or large mass ratios 
are indicated. 



host star. Apparently, for these systems, the trend observed 
by Winn et al. (2010a) still holds. There are three apparent 
exceptions to the rule: HAT-P-11, WASP-8, and HD 80606. 
These systems are special in other ways too, by virtue of hav- 
ing either an unusually low planet mass or an unusually long 
orbital period. They represent three of the four systems for 
which the orbital period is greater than 7 days or the planet 
has a mass lower than 0.2 Mj up . In this sense they least resem- 
ble the typical "hot Jupiter". We will discuss these important 
cases in the following paragraphs. 

The explanation for the relationship between r en - and A 
could fall into one of two categories: (i) the formation and 
evolution of hot Jupiters is different for hot stars than for cool 
stars, which for some reason results in higher obliquities in 
the hot stars, (ii) The distribution of obliquities is originally 
broad for both hot stars and cool stars, but they evolve differ- 
ently depending on a parameter closely associated with T e ff. 
Winn et al. (2010a) suggested that the second scenario is more 
likely and that the factor associated with temperature is the 
rate of tidal dissipation due to the tide raised by the planet. 

The reason for this suspicion was that T e ff s» 6250 K is 
not an arbitrary temperature, but rather represents an approx- 
imate boundary over which the internal structure of a main- 
sequence star changes substantially. Stars hotter than this 
level have very thin or absent convective envelopes, with the 
mass of the envelope dropping below about 0.002 M at 
6250 K (Pinsonneault et al. 2001). (For the Sun, the mass 
of the convective envelope is around 0.02 M ). 

Independently of this theoretical expectation, there is dra- 
matic empirical evidence for a transition in stellar properties 
across the 6250 K divide: hotter stars are observed to rotate 
more rapidly. In the lower panel of Figure 20, we plot the 
projected rotation speeds of a sample of ~ 1000 stars from 
the catalogue by Valenti & Fischer (2005). The projected ro- 
tational speed vsini* increases rapidly around 6250 K. For 
F0 stars, the rotation speed can approach 200 kms" 1 . It is 
thought that stellar rotation together with the convection in 



the envelope create a magnetic field coupling to the ionized 
stellar wind far beyond the stellar radius and thereby trans- 
port angular momentum away from the stellar rotation (see, 
e.g., Barnes 2003, for further discussion.) Presumably this 
magnetic braking is less effective for stars without convec- 
tive envelopes, leading to the observed rapid increase in stel- 
lar rotation speeds towards earlier spectral type. Judging from 
Figure 20, the transition from low obliquity to high obliquity 
seems to be linked empirically to this transition from slowly- 
rotating to rapidly-rotating stars. 

The presence of a convective envelope is also expected to 
change the rate of dissipation of the energy in tidal oscilla- 
tions. Energy contained in tidal bulges is thought to be more 
effectively dissipated by turbulent eddies in convective en- 
velopes than by any mechanism acting in a radiative enve- 
lope. See Zahn (2008) for a review on the theory of tidal 
interactions. Torres et al. (2010) and Mazeh (2008) review 
the evidence for tidal interactions in double-star systems, and 
provide more access points to the literature. 

If tidal evolution is responsible for the difference in stel- 
lar obliquities between cool and hot stars, then there should 
also be a correlation between the mass ratio of the star and 
planet and the degree of alignment. The higher the mass of 
the planet (m p ) relative to the mass of the star (M*), the faster 
its tides can align the stellar and orbital angular momentum 
vectors. Furthermore, one would expect an inverse correlation 
between the scaled semi-major axis (R±/a) and the obliquity. 
A more distant companion will raise smaller tides. We will 
next investigate whether such correlations exist. We will also 
investigate if the age of the systems is an important factor in 
setting the degree of alignment in these systems. 

Mass Ratio — In Figure 21 the measured projected obliqui- 
ties are plotted as a function of the mass ratio between planet 
and host star. Higher obliquities are measured for systems in 
which the mass of the planet is relatively small. This is what 
would be expected if tides are responsible for the obliquity 
distribution. Massive planets raise stronger tides. This trend 
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was observed before (e.g. Hebrard et al. 201 la), though at that 
time it was not interpreted in terms of tidal interaction. 

The interpretation is not clear from this comparison alone, 
though. Note for example that the most massive planets are 
found around hot stars, which should have weaker dissipation 
that counteracts the effect of the more massive planet to at 
least some degree. In a subsequent section we will attempt 
to take both these effects (and that of orbital distance) into 
account at once. 

Distance dependence — Does the degree of alignment depend 
on the scaled distance (a/R+)l Figure 22 gives a mixed an- 
swer to this question. Focusing on systems with cool stars 
(blue open circles) there seems to be a trend of obliquities as 
a function of the scaled distance. The data suggest good align- 
ment for all systems with a/R^ < 10. Three of the four sys- 
tems with greater distances have significant misalignments. 

No such dependence is observed for systems with stars 
which are close to 6250 K or hotter. However the range of 
distances that is spanned by the hot-star systems is very small, 
less then an order of magnitude. In contrast, the cool-star sys- 
tems probe nearly two orders of magnitude. Note that two 
of the close-in, misaligned systems are the systems with the 
hottest host stars (WASP-33 and KOI- 13). This lack of align- 
ment finds an explanation in the tidal hypothesis: despite the 
tight orbits, the tidal dissipation rate may be relatively low due 
to thin or non-existent convective layer. 

Age — Under the tidal hypothesis, older systems should tend 
to be closer to alignment than younger systems, all else be- 
ing equal. This is because in older systems, tides have had a 
longer interval over which to act. Included in "all else being 
equal" is the underlying assumption that the arrival time of 
the hot Jupiter to its close-in orbit is the same in all systems. 

Triaud (2011) presented empirical evidence that the degree 
of misalignment depends chiefly on the age of the system. 
He found that all systems in his sample with ages greater 
than 2.5 Gyr are aligned (see his Figure 2). His sample con- 
sisted only of those stars with an estimated mass greater than 
1 .2 Mq , since it is harder to determine a reliable age for lower- 
mass systems. 

Stars with a mass of £1.2 M Q develop a significant con- 
vective envelope during their main sequence lifetime, even if 
they were too hot to have a significant convective envelope 
on the "zero age" main sequence. In Figure 23 we plot the 
projected obliquities as function of stellar age for stars with 
> 1.2Mq. And indeed all the aligned and older systems 
are cool enough to have a significant convective envelope. 
Figure 23 represents, therefore, a similar pattern as seen in 
Figure 20 with a slight shift in the variable and for a subset of 
systems (only stars withM^ > 1.2M Q ). It seems as though the 
development of a convective envelope with age, rather than 
the age itself, might be driving the degree of alignment. 

5.3. Tidal time scale 

As we have seen in the last section that the degree of align- 
ment is correlated with stellar temperature, the mass ratio, and 
possibly the orbital distance. We now try in this section to 
establish a single quantitative relationship between the degree 
of alignment and those parameters. Ideally we could calculate 
a theoretical alignment timescale for each system, and com- 
pare that timescale to the estimated age of the system. We 
could then check if systems with a relatively short timescale 
(fast alignment) tend to have low obliquities, and systems with 



timescales comparable to the lifetime of the system (or larger) 
tend to have high obliquities. 

Calculating timescales needed to synchronize and align 
stellar rotation is a complex task. Apart from the parame- 
ters mentioned above, there are other parameters that would 
influence the time needed for alignment. For example the to- 
tal amount of angular momentum stored in the stellar rotation, 
and the driving frequency of the tidal force (i.e. twice the or- 
bital frequency), are expected to be important. In addition, 
the rate of dissipation is not expected to be constant over Gyr 
timescales due to the contemporaneous evolution in orbital 
distance and eccentricity, and due to stellar evolution. Even 
worse, the specific mechanisms for dissipating tidal energy 
are not completely understood, neither for stars with radiative 
envelopes nor for stars with convective envelopes. Neverthe- 
less there are some simple considerations we may employ to 
obtain approximate timescales for alignment. 

1. We can use the formulae provided by Zahn (1977) 
for synchronization. The coefficients in these formu- 
lae are difficult to derive from theory alone, but they 
can be calibrated with observations made in binary star 
systems. By observing the maximum orbital distance 
within which binary stars are observed to be spin-orbit 
synchronized, and knowing the lifetime of the stars on 
the main sequence, the relevant parameters can be esti- 
mated. To apply this to our sample two different formu- 
lae are needed. One for cool stars which have convec- 
tive envelopes (CE) and hot stars which have radiative 
envelopes (RA). Therefore this approach has the virtue 
of empirical calibration, although the calibration is for 
star-star interactions rather than planet-star interactions, 
and the calibration is for spin synchronization rather 
than reorientation. One complication is that to apply 
these formulae we have to make a binary decision on 
whether a star is "convective" or "radiative" which does 
not reflect the gradual thinning of the convective enve- 
lope with increasing stellar temperature. We choose a 
temperature of 6250 K for this boundary. 

2. Assuming that the alignment timescale due to dissipa- 
tion in convective envelopes (tqe) is always shorter than 
the time needed for alignment by forces in radiative en- 
velopes (tra) we can try to derive a simple relationship 
between the mass contained in the convective envelope 
and the alignment timescale r. This would have the ad- 
vantage that the gradual decrease of mass in the convec- 
tive envelope can be easily incorporated, but we ignore 
here any possible additional dissipation mechanism in 
the radiative envelope which for higher temperatures 
would become important. In addition it is not obvious 
why r" 1 should depend linearly on the mass contained 
in the convective envelope, nor can we be completely 
confident in our estimate of the convective mass based 
only on the observable parameters of the stellar photo- 
sphere. And of course the convective mass is not really 
a constant over Gyr timescales. 

The simplifications made by either approach should cause 
us not to expect a perfect and deterministic relationship be- 
tween our theoretical parameters and the observed obliqui- 
ties [and we direct the reader to Hansen (2012) for a different 
approach to this comparison]. For the first approach we ob- 
tain the following relationships between the system parame- 
ters and the convective and radiative timescale for alignment 
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FIG. 22. — Dependence of A on the scaled orbital distance. The same as Figure 2 1 , but now the measurements of the projected obliquity are plotted as a 
function of the semi-major axis divided by the stellar radius. While a trend with distance can be observed for planets around cooler stars, such a trend seems to be 
absent for planets around hotter stars. This might be due to the small range probed in distance by the obliquity measurements around hot stars. However because 
tidal forces are week in these stars only for the innermost massive planets (e.g. WASP-12 b and WASP-18 b) such a trend would be readily observable. Names of 
systems with scaled distances greater then 15, and for misaligned close in hot systems are indicated. 
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Here q is the planet-to-star mass ratio. For stars with convec- 
tive envelopes, synchronization is observed for binaries out 
to a scaled distance of ~ 40 during their main-sequence life, 
which we set to 10 Gyr for all 'cool' stars. For 'hot' stars, 
synchronization in binaries is observed out to ~ 6 times the 



scaled radius. Tidal damping is expected to be most efficient 
during the first quarter of their main sequence life (Zahn 1977) 
(which we set here to 5 Gyr). See e.g. Claret & Cunha (1997) 
for a more recent comparison between theory and observa- 
tions, mainly for stars of higher mass. 

In Figure 24 we show the projected obliquity versus the 
characteristic timescale needed for realignment. For this we 
used equation 2 for T e ff < 6250 K and equation 3 for T t ff > 
6250 K. (Both timescales were divided by 5 • 10 9 for normal- 
ization.) For most of the systems for which rapid alignment 
is expected, low projected obliquities are observed. For sys- 
tems where tides are expected to be too slow in aligning and 
synchronizing stellar rotation, a very broad range of projected 
obliquities (apparently random) is observed. 

We now return to the three apparent outliers from Figure 20: 
HAT-P-1 1, WASP-8, and HD 80606. While these are all stars 
with convective envelopes, the timescales for alignment are 
very long. This is because the scaled distances are greater 
than 15 (Figure 22). In the case of HAT-P-1 1 there is also 
the additional penalty from the relatively small planet mass 
(q ss 10~ 4 ). Thus, in this light, these three "outliers" are 
not exceptions; they have high obliquities because the tidal 
timescales are very long. 

There is one system which does seem to be an exception: 
HAT-P-32. The rotation axis of the star nearly lays in the 
plane of the orbit (A = 85 ± 1 .5°), while all other systems with 
a similar tidal timescale do have projected obliquities consis- 
tent with good alignment. Because of the obliquity near 90° 
tides couple only weakly to the stellar rotation (Lai 2012). 
The timescale for alignment in this system could be longer 
than estimated by our simple formula. If this is the only 
reason for the high obliquity in the HAT-P-32 system then 
we might expect to find more systems with short alignment 
timescales and similar obliquities in Figure 24, which do not 
exist in the current sample. 

The measured effective temperature in HAT-P-32 is 6207 ± 
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Fig. 24. — Measured projected obliquity as function of the alignment timescale calibrated from binary studies. The same symbols as in Figure 20 are 
used. This time the projected obliquities are shown as a function of the characteristic timescale needed to align the stellar and orbital axes. We used two different 
equations to calculate these timescales. One for stars with temperatures lower than 6250 K for which we assume that tidal dissipation happens due to Eddies 
in the convective envelope and one for hotter stars for which we assume that no convective envelope is present and alignment is due to radiative damping. The 
coefficients for these equations have been calibrated with synchronization timescales in double star systems. Note that both timescales have been divided by 
5-10 9 . 




1(T 4 10" 2 10° 10 2 10 4 10 6 

relative tidal-dissipation timescale x mcz / age 

Fig. 25. — Measured projected obliquity as a function of the alignment timescale estimated from the mass of the convective envelope divided by 
age. Similar to Figure 24. This time however a tidal timescale was calculated which depends next to the mass ratio and the scaled distance not on a calibrated 
coefficient but on the mass in the convective envelope. Further the age estimates of the systems have been taken into account for this ranking. KOI-13 is not 
shown in this plot. The mass contained in its convective envelope is small and therefore according to equation 4 it does practically not realign. 
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88 K, and therefore we used Equation 2, applicable to stars 
with convective envelopes. If we would have assumed an 
effective temperature of 6250 K, a value within the 1-er in- 
terval of the measurement, and used Equation 3 then we 
would have obtained a tra of 1.5 ■ 10 5 x 5 • 10 9 yr instead of 
tce = 1 .8 • 10 1 x 5 • 10 9 yr. This illustrates the aforementioned 
weakness of this first approach, that we have to make a binary 
decision on whether a star is "convective" or "radiative." 

Note also that KOI- 13 and XO-3 are "hot" stars which have 
significant misalignments, and yet they are found in between 
"cool" aligned systems. KOI-13 (T e ft = 8500 K) is hottest star 
in our sample and it is questionable if we can use the same 
tidal timescale for this system as for the other hot systems 
which are about 2000 K cooler (see Figure 20). 

In the second approach we build upon Equation 2. Now we 
do not use any empirical calibration. For each planet-hosting 
star we estimate the mass contained in the convective enve- 
lope, and assume that the rate of energy dissipation is propor- 
tional to this convective mass, 



t M„ \a 



(4) 



where M cz indicates the mass contained in the outer convec- 
tive zone and C is an unspecified proportionality constant with 
units g s" 1 . Our estimate for M cz is based on the measured 
r e ff. This ensures a gradual decrease of tidal forces with in- 
creasing temperature. To establish the relation between T e ff 
and convective mass we used the EZ-Web tool' 1 for stars with 
T eS < 7000 K, and the data from Pinsonneault et al. (2001) for 
hotter stars. To create Figure 25 we further divided the tidal 
timescale by the estimated main-sequence age, taking the un- 
certainty in the age estimate into account. As the ages are not 
well known, this leads to a substantial uncertainty in the po- 
sitioning of a system on the horizontal axis. On average, the 
hot stars are younger than the cool stars. Therefore the main 
effect of the division by age is a small shift of the hotter stars 
to the right side of the logarithmic plot. 

The ordering of the cool stars is not substantially changed, 
relative to Figure 24, but there are now a few hot stars 
with significant obliquities and with similar tidal timescales 
as some cool aligned systems. The biggest outlier in this re- 
spect is HAT-P-7. 

In summary, despite the shortcomings of our highly simpli- 
fied approaches to calculating the theoretical tidal timescale, 
and a few exceptional cases, we do find support for the claim 
that the obliquities in hot Jupiter systems undergo damping 
by tidal dissipation. Systems with short tidal timescales are 
predominantly well-aligned, while systems with longer tidal 
timescales display an apparently random obliquity distribu- 
tion. The implication is that the obliquities were once even 
more broadly distributed than we observe them today. Put dif- 
ferently, the "primordial" orbits of hot Jupiters (the orbits that 
existed shortly after the planets arrived close to the star) once 
had a very broad range of inclinations relative to the stellar 
equatorial plane. 

5.3.1. Angular momentum problem 

As Winn et al. (2010a) already pointed out, there is a theo- 
retical problem with invoking tides in this context. The angu- 
lar momentum in the stellar rotation compared to the angular 

6 This tool is made available by Richard Townsend under the following 
url: http : / /www .astro. wise, edu/ ~ town send 



momentum in the orbit (when the planet is close enough to 
significantly effect the stars rotation via tides) is so large that 
to synchronize and align the star the planet would surrender 
so much angular momentum that it would spiral into the star. 
For nearly all systems in our sample the orbital velocity (at 
periastron) is larger than stellar rotation velocity. This causes 
trailing tides, and angular momentum is transported from the 
orbit towards the stellar rotation, leading to decay of the or- 
bit (e.g. Levrard et al. 2009). Yet we see systems which have 
aligned axes and the planets have evidently survived. 

To address this problem Winn et al. (2010a) speculated that 
only the outer layers of the star synchronize and align with the 
orbit. In that case a smaller amount of the angular momen- 
tum would be transferred and the planetary inspiral would be 
avoided. It seems doubtful, though, that a separate rotation 
speed and rotation direction for the envelope relative to the 
stellar interior could be maintained for billions of years. 

More recently Lai (2012) suggested that the angular mo- 
mentum problem is not as serious as it might seem. Given 
the complexities of tidal dynamics, he argued that there is no 
strong theoretical reason why the timescale for realignment 
must equal the timescale for synchronization, and indeed he 
provided a particular theoretical tidal model in which those 
timescales can differ by orders of magnitude. In his scenario 
the planets would first align the stellar rotation, and only much 
later speed up the rotation and spiral inward. 

In this respect it is interesting that the tidal timescale cali- 
brated via synchronization apparently sorts the systems con- 
sistently relative to each other, but the overall timescale is too 
long by orders of magnitude. As mentioned above we di- 
vided the timescales displayed in Figure 24 by 5 • 10 9 . This 
could imply that realignment happens on a shorter timescale 
than synchronization. However the calibration of the synchro- 
nization timescale was done with binary star systems having 
q R3 1, and the tidal mechanism itself might be different for 
different regimes of its strength (Weinberg et al. 2012). 

One might be able to test the hypothesis of Lai (2012) by 
seeking evidence for excess rotation in stars that are thought 
to have been realigned (A ss 0°). This could be done by mea- 
suring the stellar rotation period or vsini* (if one is willing to 
assume sin/* is near unity in such systems). If an age estimate 
is also available, then one could employ the same approach as 
Schlaufman (2010) to assess whether or not the star is rotating 
at a typical rate, or if it is in the process of being spun up by 
the planet. This type of analysis was pursued by Pont (2009), 
though not with this specific hypothesis test in mind. This 
analysis could be profitably revisited now that many more 
systems are available for study. There is one caveat to this ap- 
proach, which is that stars undergo very rapid spin evolution 
early in their lives due to disk-locking and magnetic braking, 
i.e., for reasons unrelated to planets. If hot Jupiters arrive very 
early in the star's history, the realignment might happen in an 
epoch of rapid decrease of stellar rotation and any memory of 
an increased rotation due to tides raised by the planet might 
be lost. 

The recent work by Hansen (2012) presented a calibration 
of the equilibrium tide theory using the measured parameters 
of hot- Jupiter systems. While the approach taken in his work 
is different from ours, he arrived at similar conclusions to 
those described here: tides have shaped the obliquity distri- 
bution in these systems, and there is no theoretical need for 
core-envelope decoupling. 
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5.4. High obliquities: a result of dynamical interactions, or 
initially inclined disks? 

To interpret the finding that the host stars of hot Jupiters 
once had a broad distribution of obliquities, we must answer 
a crucial question. We need to know if the original obliquity is 
related to the existence of the hot Jupiter, or if stars and their 
protoplanetary disks are frequently misaligned for reasons un- 
related to hot Jupiters. One might expect an initially close 
alignment between a star and its protoplanetary disk, as is ob- 
served in the Solar system and has been generally assumed in 
the exoplanet literature. However this is not a foregone con- 
clusion, and indeed several authors have recently challenged 
this assumption, proposing that the Sun's low obliquity may 
be an atypical case. 

Bate et al. (2010) proposed that a disk might become in- 
clined with respect to the rotation axis of the central star, as 
a result of the complex accretion environment within a star 
cluster. In such a dense environment the tidal interaction with 
a companion star or other nearby stars could produce chaotic 
perturbations in the orbits of infalling material, with the mate- 
rial accreting later (destined to become planets) having a dif- 
ferent orientation than the material that accreted earlier onto 
the star. Thies et al. (201 1) studied the process of inclined in- 
fall of gas in detail and found that short period planets on ec- 
centric and inclined orbits can be created in this way. A com- 
pletely different mechanism for generating primordial mis- 
alignments was proposed by Lai et al. (2011), relying on a 
magnetic interaction between a young star and the inner edge 
of its accretion disk. 

In these scenarios, the star has a high obliquity even though 
the planets may have never left the plane of the disk out of 
which they have formed, and therefore the measurements of 
obliquities bear information about the processes surrounding 
star formation rather than planet migration. How can one dis- 
tinguish between misalignment created during the time the 
planet is still embedded in the disk or after the disk dissipated? 

One approach, pursued by Watson et al. (2011), is to as- 
sess the degree of alignment between stars and their debris 
disks. Assuming that the stars as well as the debris disks trace 
the alignment of their predecessors one would learn about the 
degree of the alignment during the final stages of planet for- 
mation. Watson et al. (2011) found the inclinations of debris 
disks and their stars in a sample of 9 systems to be consis- 
tent with good alignment. They caution that in their study 
only systems with r e g < 6140 K have been observed and that 
misaligned system are found around hotter stars. However as 
we have argued above, the found low obliquities that prevail 
around cool stars may be a consequence of tidal evolution and 
not of the mechanism which creates the obliquities. 

Another approach is to measure the obliquities in binary 
star systems. If disks would be tilted relative to stellar rota- 
tion due to close encounters, then this could also lead to tilted 
rotation axes in double star systems. There should also be 
a trend towards misalignment with larger separation between 
the components in these systems. Conducting such measure- 
ments and seeking evidence for such trends is one of the goals 
of an ongoing observational program entitled BANANA (Bi- 
naries Are Not Always Neatly Aligned; Albrecht et al. 201 la). 

A more direct way to answer the question raised in the pre- 
ceding section has recently become possible, thanks to the dis- 
covery of systems with multiple transiting planets. A number 
of arguments support the idea that the orbital planes in such 
systems are closely aligned; most recently Fabrycky et al. 



(2012) used the measured transit durations to show that the 
typical mutual inclinations are of order 2°. Therefore it is rea- 
sonable to assume that the multiple planetary orbits trace the 
plane in which the planets formed. Any disruptive dynamical 
interactions, such as those which have been proposed to ex- 
plain hot Jupiters, would likely have produced higher mutual 
inclinations in the multiple-transiting systems. Under that as- 
sumption, RM measurements (or other measures of obliquity) 
in those multiple-transiting systems would establish the angle 
between the circumstellar disk and the stellar equator. If good 
alignment is found to be the rule, then the high obliquities 
in hot Jupiter systems would be more readily interpreted as a 
consequence of planet migration than as primordial star-disk 
misalignments. 

6. SUMMARY 

In this paper we have presented new observations of the RM 
effect for 14 systems harboring hot Jupiters. In addition we 
critically reviewed the literature, in some cases re-analyzing 
data that had been obtained previously in order to conform 
with our protocols. We then used these data to show that the 
distribution in obliquities is consistent with being shaped by 
tides raised by the hot Jupiters on the stars. For this we revis- 
ited the correlation between the projected obliquity and the 
effective temperature discovered by Winn et al. (2010a), now 
with a sample of RM measurements twice as large as was then 
available. We showed that the new measurements agree with 
the pattern proposed by Winn et al. (2010a). With the en- 
larged sample we showed that obliquity in systems with close 
in massive planets further depend on the mass ratio and the 
distance between star and planet, in roughly the manner ex- 
pected if tides are responsible for the low obliquities. 

Motivated by these results we then devised two different 
parameters that represent, at least crudely, the theoretical tidal 
timescales. This showed that systems which are expected to 
align fast are all showing projections of the obliquities which 
are consistent with good alignment. In contrast, systems for 
which tidal interaction is expected to be weak, due to the stel- 
lar structure, distance, or mass ratio, show a nearly random 
distribution of projected obliquities. Our interpretation is that 
stars with hot Jupiters once had a very broad range of obliq- 
uities. It is tempting to argue further that the large obliquities 
originate from the same process that produces hot Jupiters, 
thereby favoring explanations involving dynamical scattering 
or the Kozai effect, and disfavoring the gradual inspiral due 
to torques in a protoplanetary disk. However, more observa- 
tions are needed to check on the possibility that stars and their 
disks are frequently misaligned for reasons unrelated to hot 
Jupiters. Among these observations are the extension of RM 
studies to planets other than hot Jupiters, and measurements 
of obliquities in binary star systems and in systems with mul- 
tiple transiting planets. 
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